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<N Abstract 



The generic real {(3 = 1) and complex {(3 = 2) two-qubit states are 9-dimensional and 15- 



O 

dimensional in nature, respectively. The total volumes of the spaces they occupy with respect to 



the Hilbert-Schmidt and Bures metrics are obtainable as special cases of formulas of Zyczkowski 
and Sommers. We claim that if one could determine certain metric- independent 3-dimensional 
"ei^enva/Me-parameterized separability functions" (EPSFs), S^'^\\i . . . A4), then these formulas 
qh could be readily modified so as to yield the Hilbert-Schmidt and Bures volumes occupied by only the 

separable two-qubit states (and hence associated separability probabilities). Motivated by analogous 

m 

earlier analyses of ll diagonal-entry-pa,ia,meterized separability functions", we further explore the 
possibility that such 3-dimensional EPSFs might, in turn, be expressible as univariate functions of 
some special relevant variable-which we hypothesize to be the maximal concurrence (0 < C < 1) 
over spectral orbits. Extensive numerical results we obtain are rather closely supportive of this 

o 



hypothesis. Both the real and complex estimated EPSFs exhibit clearly pronounced jumps of 
magnitude roug hly 50% at C = |, as well as a number of additional matching discontinuities. 
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I. INTRODUCTION 



In a pair of major, skillful papers, using concepts of random matrix theory, Zyczkowski 
and Sommers were able to obtain exact formulas for the total volumes-both in terms of the 
Hilbert-Schmidt (HS) metric [1] and Bures (minimal monotone) metric [2] -of the (N 2 — 1)- 
dimensional convex set of N x N complex density matrices and the ({N 2 + N — 2)/2)- 
dimensional convex set of N x N real density matrices, representing iV-level quantum sys- 
tems (cf. [3] [U sees. 14.3, 14.4]). In two recent studies, we have been interested in the 
question of how to modify/truncate, in some natural manner (by multiplying certain in- 
tegrands by relevant functions), these formulas of Zyczkowski and Sommers, so that they 
will yield not total volumes, but only the (lesser, strictly included) volumes occupied by the 
separable /nonentangled states [HIE] (cf. [?]). We will below report some interesting progress 
in this regard, in relation to the two-qubit (N = 4) states. 

To begin, we present two parallel formulas from [I] and [2] for certain generalized nor- 
malization constants (Cjy ) used in the total HS and Bures volume computations. (Some 
notation and formatting has been altered.) For the HS case, we have [U eq. (4.1)] (cf. jH 
eq. (14.35)]), 
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and for the Bures case eq. (3.19)] (cf. (3 eq. (14.46)]), 
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The A's are the N (nonnegative) eigenvalues-constrained to sum to 1-of the corresponding 
NxN density matrices, while the parameter (3 is a "Dyson index", with (3=1 corresponding 
to the real case, and (3 = 2, the complex case (and (3 = 4, the quaternionic case, not explicitly 
discussed in [HE]). The parameter a will be equal to 1 for the case-of immediate interest 
to us here-of generically nondegenerate density matrices. 

A. Objective 

Our goal, in overall terms, is to find metric-independent (separability) functions, 

(Ai . . . Ajv) , (3) 
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which, if inserted into formulas ([TJ and ^ under the integral signs, as simple multiplicative 
factors, will yield separable-rather than total-volumes when the resulting modified C^f s 
are employed in exactly the same auxiliary computations (involving flag manifolds) in pQ 
and [2j as the C^'s given by ([TJ and ^ were there. More specifically here, our numerical 
analyses will be restricted to the N = 4 and (3 = 2 (complex) and (3 = 1 (real) cases. 

Our metric-independent goal is plausible for the following reason. Precisely the same 
preliminary integrations-respecting the separability constraints-over the non-eigenvalue pa- 
rameters [possibly, Euler angles [E] [13 App. I]] must be performed for both metrics before 
arriving at the stage at which we must concern ourselves with the remaining integration 
over the eigenvalues and the differences that are now clearly apparent between metrics in 
their corresponding measures over the simplex of eigenvalues. Although we are not able to 
explicitly/symbolically determine what the results of these preliminary integrations might 
be (the computational challenges are certainly considerable), they must-whatever form 
they may take-obviously be the same for both metrics in question. Our goal here is to 
understand-with the assistance of numerical methods-what functional forms these prelimi- 
nary (12-dimensional in the complex case, and 9-dimensional in the real case) integrations 
yield. 

B. Maximal concurrence and absolute separability 

The further narrower specific focus of this study will be to explore the possibility that 
there exists a functional relationship of the form, 

si m (X 1 ...X A ) = a^\C(X 1 ...X A )), (4) 

where a^\x) are some unknown univariate (one-dimensional) functions and 

C(Ai . . . A 4 ) = max{0, Ai - A 3 - 2 V%M, Ai > A 2 > A 3 > A 4 , (5) 

is the maximal concurrence over spectral orbits of two-qubit density matrices [9j sec. VII] 

maun. 

For two-qubit states, C G [0,1], with C = corresponding to the absolutely separable 
states. That is, no density matrix with C = can be nonseparable/entangled [12] • (In a 
recent study, we were able to obtain exact expressions-involving the tetrahedral dihedral 
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angle cos -1 (|)-for the contributions to the Hilbert-Schmidt real and complex two-qubit 
volumes for those states with C = 0, and to numerically estimate the Bures counterparts 
[HI sees. III.B, III.C]. In numerical terms, the HS absolute separability probability of generic 
complex two-qubit states is 0.00365826, and the Bures counterpart, 0.000161792. The HS 
real analogue is 0.0348338.) The concurrence itself is a widely used entanglement measure 
of bipartite mixed states jH eq. (15.26)]. 



C. Motivation 

Certainly part of our motivation for advancing the ansatz Q was that an analogous mod- 
eling of a trivariate function in terms of a univariate function was found to hold-making use 
of the Bloore (correlation coefficient) parameterization of density matrices [13]-for diagonal- 
enin/-parameterized separability functions P~H eq. (6)] [T5]. This led to substantial insights- 
and exact conjectures and ^)-with regard to Hilbert-Schmidt (complex and real) two- 
qubit separability probabilities. (The Dyson indices (3 played a central analytical role there, 
in relating real and complex [and quaternionic] results, but not apparently-as far as we can 
perceive-in the analyses to be presented below.) 



II. NUMERICS 



A. Methodology 



We do find encouragement in advancing the ansatz Q by the extensive numerical results 
we generate, in that our estimates of crW(C) and a^^C) shown in Fig. [I] rather closely 



reproduce-as we will indicate below (sec. |II B (-other (independent) numerical results and 



accompanying conjectures we have previously obtained. 

The (3 = 2 complex curve in Fig. [T] shown is based on the use for quasi-Monte Carlo 
numerical integration of 26,300,000 12- dimensional [Tezuka-Faure (TF) [16]] low- discrepancy 
points, and the (3 = 1 case, on 33,000,000 6-dimensional such points. (The TF procedure- 
programmed in Mathematica by Giray Okten [17] -is not conducive to the placing of error 
bars on the results, though later routines developed by him are.) These points comprise 
sample values, respectively, of the 12 Euler angles used to parameterize S77(4) and the 6 
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max. cone. 



FIG. 1: Joint plot of estimated (real [blue, = 1] and complex [red, (3 = 2]) functions of the 
maximal concurrence over spectral orbits, Sf'^Ai . . . A 4 ) = a^\C{\ . . . A 4 )). Note evident jumps 
in both functions when the maximal concurrence equals 0.5. The graphs cross at C = 0.181245. 
For C = 0, <j(C) = 1, so all associated density matrices are separable. 



Euler angles used for 5*0(4). For each TF-point, 499 auxiliary computations were carried 
out-in addition to that of the corresponding Haar measure associated with the Euler angles- 
for sets of eigenvalues with values of maximal concurrence running at equally-spaced intervals 
from ± to fg . 

Each density matrix generated-corresponding to a specific set of eigenvalues with fixed C 
and Euler angles [8] [6j App. I]-was tested for separability. Prior to the quasi-Monte Carlo 
runs, we established a database-using the Mathematica command "FindInstance"-of 100 
sets of four eigenvalues for each of the equally-spaced 499 values of C . One of the 100 sets 
was randomly selected [and then randomly permuted] for each of the TF-points and each of 
the 499 iterations. This "random generation" of sets of eigenvalues with fixed values of C 
is clearly less than an ideal procedure, but it was what we found to be practical under the 



particular circumstances. (In sec. HE, we manage to improve upon this approach.) 



Several weeks of MacMini computer time were used for each of the two sets-real and 
complex-of calculations. (Along with the computations concerning the maximal concurrence 
(5 ), we also carried out a fully parallel set of computations using the related variable, 2x " X ~' X| 



Ai-a 3 



Those results, however, seemed comparatively disappointing in their predictive power, so we 
do not detail them here.) 



5 



B. Evaluation of numerical results 

Let us now appraise our estimated functions (Fig. [T]) by seeing how well they are able 
to reproduce previous related results, themselves based on very extensive analyses (mostly 
involving quasi-Monte Carlo integration also). 

1. Complex case 

Use of the complex {(3 = 2) function in Fig. [I] impressively explains 98.7253% of the 
variance of the estimated trivariate eigenvalue-parameterized separability function for C > 
presented in [5j sec. III.B]. It also yields an estimate of 0.254756 for the Hilbert- Schmidt 
separability probability, while our exact conjecture from [15] is ^ ~ 0.242424. Further, the 
Bures separability probability estimate yielded is 0.0692753, while our conjectured value is 
1M °(; g 1+ ^) « 0.Q733389 gB|. 

2. Real case 

Passing to the real (/3 = 1) case, we had previously formulated the conjecture that the 
HS separability probability is A ~ 0.470588 [15], sec. 9.1]. Our estimate based on the (blue) 
function shown in Fig. [l]is 0.480302. (The corresponding estimate-for which we have no prior 
conjecture-for the Bures real two-qubit separability probability is 0.212152.) Further, we 
are able to reproduce 97.7502% of the variation in the corresponding trivariate function for 
C > 0. (This last function had been estimated using a recent Euler-angle parameterization 
of 5*0(4), obtained by S. Cacciatori [6, App. I]. It was derived by Cacciatori after the 
submission of [5], and thus not reported nor used there, although its complex counterpart- 
based on 3,600,000 Tezuka-Faure points-had been [5j sec. III.B], while the real case was 
based on a considerably lesser number of TF-points, 700,000.) 

C. Jumps near C = \ 

For the real (/? = 1) case, the jump near C(Ai . . . A4) = | is from approximately 0.118696 
to 0.180357, and in the complex (/3 = 2) case, from 0.0439255 to 0.651586. The magnitudes 
of the two jumps are then quite comparable, being respectively, 51.964% and 51.488%. In 
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0.510 



FIG. 2: Same plot as Fig. [IJ restricted to the vicinity of C = |, and the complex (red) curve being 
amplified by a factor of 2.8 



o-'(C) 




max. cone. 



FIG. 3: Joint plot of derivatives with respect to C of the estimated (real [blue, (3=1] and complex 
[red, (3 = 2]) functions of the maximal concurrence over spectral orbits, a^\C(X . . . A4)). Spikes 
are observable-froi/i for the real and complex cases-at: \ = 0.5; ItK = 0.294; Mk. = 0.204; and 
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0.34. 



Fig. [2J we replot the curves shown in Fig. [T] in the immediate vicinity of C = | , but amplify 
the complex (red) curve by a factor of 2.8. We perceive a very close similarity in shape. 



D. Additional discontinuities 

In Fig. [3] we show the derivatives with respect to C of the estimates of cr^'(C). (Fig. [4] 
is a plot of the same two curves, except that we have now added 10 to the derivative in the 
real case and subtracted 10 in the complex case, so that the discontinuities can be more 
readily distinguished and compared.) In addition to the already-discussed behavior at 
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FIG. 4: Same as Fig. [3] except that the real (blue) curve has been translated upwards by 10 units 
and the complex (red) curve downwards by 10 units, so that the individual discontinuities in the 
two derivatives can be more readily seen and compared. 

C = 1/2, we see-both in the real and complex cases-a secondary spike at |^ = 0.294, and 
lesser spikes at ^ = 0.204 and |^ = ^ = 0.34. So, all the observed spikes, signaling 
what we presume are discontinuities in the cr^'s, and concomitant nontrivial piecewise 
behavior-indicative of different separability constraints becoming active/binding or not-are 
for C < \. The point C = §^ = 0.102 may also be a discontinuity, at least in the complex 
case. We could detect no apparent spikes/discontinuities in the upper half- range, C G [|, 1]. 
(In a somewhat analogous study of two-qubit i/jree-parameter HS separability probabilities, 
intricate piecewise continuous behavior [involving the golden ratio] was observed [TjJl eq. 
(37) and Fig. 11].) 

In Fig. [5] we show the segments of the estimated functions cr^(C) between the two 
discontinuities, ^ = 0.204 and |^ = 0.34. The behavior seems very close to linear for both 
curves, except for the intermediate discontinuity at ^ = 0.294. 



E. Supplementary analyses 

Since the completion of the extensive numerical analyses described above, we have un- 
dertaken supplementary, parallel analyses in which 5,000 (rather than 500) sub intervals of 
C G [0, 1] are employed, as well as an improved method is used for sampling random eigen- 
values with fixed values of C (using the Mathematica Findlnstance command, now with a 
random seed). These results so far seem largely consistent with those already described. 
However, the new analogues of the plots of derivatives, Figs. [3] and |4| are still much too 
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FIG. 5: Joint plot of estimated (real [blue, (3 = 1} and complex [red, (3 = 2]) functions of the 
maximal concurrence over spectral orbits, S^'^\\\ . . . A4) = a^\C{\ . . . A4)). The graphs are 
obviously close to linear between the discontinuities = 0.204 and ^ = 0.34, except for the 
intermediate discontinuity at ^ = 0.294. To high accuracy, the real (blue) curve can be fitted by 
the line 1.07614 - 1.99472C and the complex (red) curve, by 1.19822 - 2.69548C. 



cr(C) 




FIG. 6: The same plot as Fig. [2] but based on our ongoing supplementary analysis with finer 
resolution in C and enhanced eigenvalue sampling. The twin jumps in the estimated eigenvalue- 
parameterized real and complex separability functions near C = \ are now certainly indisputably 
clear. 

rough in character to detect the presence of any secondary (non-jump) discontinuities. But, 
even at this stage (having tested 30,400 times the separability of 4,999 complex density 
matrices, and 28,100 times the separability of 4,999 real density matrices), we can produce 
the interesting counterpart (Fig. |6j) to Fig. [2J in which the two jumps near C = | of roughly 
50% magnitude are clearly unmistakable. Further, the highly linear behavior displayed in 
Fig. [5] also reappears (Fig. [7]). 
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FIG. 7: Same plot as Fig. [5] but based on our ongoing supplementary analysis with finer resolution 
in C and enhanced eigenvalue sampling 

III. CONCLUDING REMARKS 

For the real and complex two-qubit systems, we have investigated the possibility that the 
associated (t/iree-dimensional) ezgenfa/t/e-parameterized separability functions-conceptually 
important for computing separability probabilities— axe expressible as one-dimensional func- 
tions (o"(C)) of the maximal concurrence over spectral orbits (C G [0, 1]). Our numerical es- 
timates, in this regard, have been encouraging, in that they closely reproduce independently- 
generated numerical results and exact conjectures concerning separability probabilities based 
on the Hilbert-Schmidt and Bures (minimal monotone) metrics over the two-qubit systems, 
and based on the use of diagonal- en£n/-parameterized separability functions. Plots of the 
real and complex versions of c(C) both exhibit jumps of approximately 50% magnitude 
near the midpoint, C = |, and both also indicate the presence of, at least, three further 
(non-jump) discontinuities (C ~ 0.204, 0.294, 0.34), apparently indicative of points at which 
certain distinct separability constraints become either active/binding or not. Over the in- 
terval C G [0.204, 0.34], the real and complex fitted functions c(C) both appear to be simply 
linear (except at C ps 0.294). 

We have principally studied above the possibility (4) that the ostensibly trivariate two- 
qubit eigenvalue-separability functions can be equivalently expressed as univariate functions 
of only a single variable, that is, the maximal concurrence C over spectral orbits [9] . Since we 
have unfortunately not been able to fully formally resolve this issue-although our supporting 
evidence for this proposition is intriguing-we can not also presently fully eliminate the possi- 
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bility that one or even two (yet unspecified) variables supplemental to C are, in fact, needed, 
and that the corresponding separability function is not in fact strictly univariate in nature 
(as we do know it definitely is the case with the two-qubit diagonal-entry-parameterized 
separability functions [15]). 

It presently appears somewhat problematical to extend the line of analysis above to 
the qnbit- qutrit (N = 6) case. In addition, to simply the greatly increased computational 
burden that would be involved, there does not seem to be a maximal concurrence formula 
comparable to the two-qubit one ^ with the requisite properties we have utilized [HI p. 
102108-16]. 

We have examined the relationship between separability and entanglement in a specific 
analytical setting-involving eigenvalue-parameterized separability functions and the use of 
the Zyczkowski-Sommers formulas pfl |2] for the total Hilbert-Schmidt and Bures volumes 
of the N x N density matrices. A number of studies of Batle, Casas, A. Plastino and A. 
R. Plastino (for example, [20]) have also focused on the relationship between separability 
and entanglement, but in somewhat different analytical frameworks (typically involving the 
ZHSL measure [7], which is uniform over the eigenvalue simplex). The closest we can come, it 
seems, to a direct comparison with their analyses, is to note that the dot-dashed curve in Fig. 
2 of [20] is based on the Hilbert-Schmidt metric, and their x-axis is the Bures distance, while 
we have employed the maximal concurrence C on the x-axis in the somewhat comparable 
Fig. 1 above (cf. [21]). Both theirs and our plots are, in general, downward-decreasing, 
but theirs gives no indication of any discontinuities. Also, their plot is of the separability 
probability, while ours is of the (presumed univariate) eigenvalue-parameterized separability 
function, to be used in the computation of the probability. 
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